1 An Overview of Robotic Mechanical Systems 


1.1 Machine: A historical account 


Here is an account of the definitions of machine, taken from (Dudita and Diaconescu, 1987): 


Different definitions of machine have been given by scholars for more than two millennia, starting with 
Vitruvius in 28 B.C., namely, 


e A machine is a combination (system, assemblage) of moving material bodies (Vitruvius, 28 B.C.; 
Hachette, 1811; Borgnis, 1818; Beck, 1859; Reuleaux, 1875; Koenigs, 1901) 


e A machine is generally composed of three parts: a motor part, a transmission part, and an exe- 
cution part (Euler, 1753; Bogolyubov, 1976) 


e A machine produces mechanical work, or performs productive operations, actions, or effects (Vit- 
ruvius, 28 B.C.; Poncelet, 1824; Reuleaux, 1900; Koenigs, 1901; Bogolyubov, 1976) 


e A machine transforms or transmits forces (Vitruvius, 28 B.C.; Leupold, 1724; Euler, 1753; Bo- 
golyubov, 1976; Reuleaux, 1900; Koenigs, 1901) 


e A machine is characterized by deterministic motions (Hachette, 1811; Leupold, 1724; Reuleaux, 
1875; Borgnis, 1818; Reuleaux, 1900) 


e A machine is an artifact (Leupold, 1724) 
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1.3 Definitions for “machine,” “mechanism,” and “linkage”: 
i ? 


Machine 


e Definitions in Merriam Webster’s Collegiate Dictionary (on-line, 2002): 
— (archaic): a constructed thing whether material or immaterial. 
— an assemblage of parts that transmit forces, motion, and energy one to another in a prede- 
termined manner 
— an instrument (as a lever) designed to transmit or modify the application of power, force, or 
motion 
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— a mechanically, electrically, or electronically operated device for performing a task (a calcu- 
lating machine, a card-sorting machine) 


Comment: comprehensive definitions when considered as a whole 

e An apparatus for transformation of power, materials, and information to substitute or simplify 
physical or intellectual work (Frolov, 1987). Comment: a comprehensive definition, that includes 
computers 

e Mechanical system that performs a specific task, such as the forming of material, and the trans- 


ference and transformation of motion and force, Vol. 38, Nos. 7-10 (2003) of Mechanism and 
Machine Theory on Standardization of Terminology. Comment: leaves computers out 


e An apparatus for applying mechanical power, having several parts, each with definite function 
(The Concise Oxford Dictionary). Comment: same as above 


e An apparatus consisting of interrelated parts with separate functions, used in the performance of 
some kind of work (The Random House College Dictionary). Comment: ditto 


e Any system in which a specific correspondence exists between an input form of energy or in- 
formation and the corresponding ones at the output (Loosely translated from Le Petit Robert). 
Comment: as comprehensive as Frolov’s 


Mechanism 


e A piece of machinery (Merriam Webster’s Collegiate Dictionary (on-line, 2002)). Comment: too 
vague 

e Definitions in Vol. 38, Nos. 7-10 (2003) of Mechanism and Machine Theory on Standardization 
of Terminology. 

— System of bodies designed to convert motions of, and forces on, one or several bodies into 
constrained motions of, and forces on, other bodies. Comment: English could be terser, but 
idea is fine. 

— Kinematic chain with one of its components (link or joint) connected to the frame. Comment: 
confuses mechanism with its representation as a kinematic chain 

e Structure, adaptation of parts of machine; system of mutually adapted parts working together 
(as) in machine (The Concise Oxford Dictionary) 

e An assembly of moving parts performing a complete functional motion (The Random House 
College Dictionary) 


e A combination layout of pieces or elements, assembled with the goal of (producing) an operation 
as a unit (Loosely translated from Le Petit Robert) 
Comment: In all above definitions, the concept of goal or task is present 


Linkage 


e Definitions in Merriam Webster’s Collegiate Dictionary (on-line, 2002): 


— a system of links. Comment: concise and comprehensive 


— a system of links or bars which are jointed together and more or less constrained by having a 
link or links fixed and by means of which straight or nearly straight lines or other point paths 
may be traced. Comment: unnecessarily cumbersome and limited to path-generating linkages 


e Kinematic chain whose joints are equivalent to lower pairs only (Vol. 38, Nos. 7-10 (2003) of 
Mechanism and Machine Theory on Standardization of Terminology). Comment: confuses linkage 
with its representation 


The Concise Ozford Dictionary of Current English, 1995, Clarendon Press, Oxford 
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Frolov, K. V., (editor), 1987, Teoriya Mechanismov i Mashin (Theory of Mechanisms and Machines), 
Vyschaya Shkola, Moscow (in Russian) 


Mechanism and Machine Theory on Standardization of Terminology, 2003, Vol. 38, Nos. 7—10 
Le Petit Robert, 1994, Dictionnaires Le Robert, Paris 

Random House Webster’s College Dictionary, 1997, Random House, New York 

Merriam Webster’s Collegiate Dictionary, on-line 2002 


1.9 Here we want to estimate the time required to multiply two floating-point numbers. Since this time is 
very small, it is more suitable to perform a large number of multiplications, say 107, which will require 
a total time of the order of seconds. 


e The Microsoft Visual C++6.0 program below was run on a Pentium IV 2.0: 


#include <iostream> 

#include <time.h> 

int main() 

{ 
clock_t start, end; 
start =clock(); 
float a, i, CLOCKS_PER_mSEC=CLOCKS_PER_SEC/1000; 
for(i = 1; i <= 10000000; i++) 
a=5 * 5; 
std::cout << ‘‘a= ‘* << a << std::endl; 
end = clock(); 
long duration=(long) (end-start) /CLOCKS_PER_mSEC; 
std::cout << ‘‘Time: ‘^ << duration<<’’ms’’<<std::endl; 
return 0; 


} 


A Pentium IV 2.0 processor required 125 ms to perform the 107 multiplications, which amounts 
to 1.25 x 1078 s/mult, or 8 x 107 mult/s. 


e The C program below was run for the same purpose on the CLUMEQ (Consortium Laval UQAM 
McGill and Eastern Quebec for High Performance Computing) supercomputer, AMD Athlon 
1900+ cluster: 


include "stdio.h" 
include "math.h" 

include <sys/time.h> 
include <sys/resource.h> 
# define RUSAGE_SELF 0 

\* calling process *\ 
main() 


{ 


# 
# 
# 
# 


int getrusage(int who, struct rusage *rusage) ; 
long diffsec, diffmsec; 

float a; 

int i; 

struct rusage time_begin, time_end; 

struct rusage *pb, *pe; 

pb=&time_begin; 
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getrusage (RUSAGE_SELF, pb); 
for(i=1; i<=100000000; i++) 
{ 
a = 5*5; 
} 
printf ("a=Zf \ 
t\n",a); 
pe=ktime_end; 
getrusage (RUSAGE_SELF , pe) ; 


diffsec=time_end.ru_utime.tv_sec-time_begin.ru_utime.tv_sec; 
diffmsec=(time_end.ru_utime.tv_usec-time_begin.ru_utime.tv_usec) ; 
printf ("User time used:%4d microsec\n",diffmsec) ; 


The CLUMEQ supercomputer required 40 ms to perform the 10” multiplications, which amounts 
to 4 x 107°? s/mult, or 2.5 x 108 mult/s. 


e On June 20, 2011, Fox News reported that Kei—Japanese for 10!®—was capable of 8 x 1015 
flops'—or eight petaflops—per second. 


1.10 Here we want to estimate the time required to add two floating-point numbers. Since this time is very 
small, it is more suitable to perform a bunch of additions, say 107, which should amount to a total 
time of the order of seconds. 


e The Microsoft Visual C++6.0 program below was run on a Pentium IV 2.0: 


#include <iostream> 
#include <time.h> 
int main() 


{ 
clock_t start, end; 
start = clock(); 
float a, i, CLOCKS_PER_mSEC=CLOCKS_PER_SEC/1000; 
for(i=1; i<=10000000; i++) 
a = 5+5; 
std::cout << "a= "<< a << std::endl; 
end = clock(); 
long duration = (long) (end-start) /CLOCKS_PER_mSEC ; 
std::cout << "Time: "<< duration<<"ms"<<std::endl1; 
return 0; 
} 


A Pentium IV 2.0 processor required 125 ms to perform the 10’ additions, which amounts to 
1.25 x 1078 s/add, or 8 x 10’ add/s. 


e The C program below was run for the same purpose on the CLUMEQ supercomputer, AMD 
Athlon 1900+ cluster: 


#include "stdio.h" 
#include "math.h" 


1See the Index for a definition of flop. 
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#include <sys/time.h> 
#include <sys/resource.h> 
#define RUSAGE_SELF 0 /* calling process */ 
main() 
{ 
int getrusage(int who, struct rusage *rusage) ; 
long diffsec, diffmsec; 
float a; 
int i; 
struct rusage time_begin, time_end; 
struct rusage *pb, *pe; 
pb=ktime_begin; 
getrusage (RUSAGE_SELF, pb); 
for(i=1; i<=10000000; i++) 
{ 
a = 5+5; 
} 


printf ("a=%f \t\n",a); 
pe=&time_end; 
getrusage (RUSAGE_SELF , pe) ; 


diffsec=time_end.ru_utime.tv_sec-time_begin.ru_utime.tv_sec; 
diffmsec=time_end.ru_utime.tv_usec-time_begin.ru_utime.tv_usec; 
printf ("User time used:%d microsec\n",diffmsec) ; 


The CLUMEQ supercomputer required 40 ms to perform the 107 additions, which amounts to 
4 x 107° add/mult, or 2.5 x 10° add/s. 


1.11 What we must find here is the largest floating-point number e that, when added to a given number a, 
leaves this number unchanged, i.e., 
ate=a 


e The Microsoft Visual C++6.0 program below was run on a Pentium IV 2.0 processor at 550 MHz: 


#include <iostream> 
//#include <time.h> 
int main() 
{ 
long i; 
//double a=1; 
float a=1; 
for(i=1; i<=10000000; i++) 
{ 
a = a/10; 
if (a==0) 
break; 
} 
std::cout << "smallest floating-point"<< i-1 << std::endl; 
return 0; 
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The value of e reported in single-precision arithmetic was 1.0e — 045; in double-precision...... 
1.0e — 323. 


e The program below was run on the CLUMEQ supercomputer, AMD Athlon 1900+ cluster: 


#include "stdio.h" 
#include "math.h" 
#include <sys/resource.h> 
main() 
{ 
int i; 
float a=1; 
for(i=1; i<=10000000; i++) 
{ 
a=a/10; 
if (a==0) 
break; 


} 


printf ("smallest floating-point=%d \n",i-1); 


This machine reported exactly the same values as its Pentium counterpart. 


2 Mathematical Background 
2.6 Let w be perpendicular to v (wTv = v’ w = 0). Since we can find two such vectors, let us call them 
w, & w2 and, moreover, let us define w; L wa, i.e., wi wo = 0. 


We have 
T =1 + uv? 


Multiplying both sides by w;, we obtain 
Tw; = w; +uv' w; = wi, ¿= 1,2 
0 


Therefore, w, & wə are eigenvectors of T, of eigenvalues à = àz = 1. 


Moreover, let 


u 
Wa = — 
llull 
Then, 
Tw3 = w3 + uv! ws = Tal + wT = Tal + vu) => Tw;= a ce vľu)ws 
u u u 


and ws is the third eigenvector of T, of eigenvalue 43 = 1 + vfu. 
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2.7 The determinant of a matrix is equal to the product of its eigenvalues. Thus 
det(1 + uv”) = det(T) = \1A2A3 =1+v7u=1+u-v 


2.8 From Q = RıR2, we have 
QQ’ = Ri R2(Rı R2)" = R; RoR? RT = RR? =1 
—— 
1 


and thus Q is orthogonal. Moreover, we have 
det(Q) = det(Rı)det(R2) = (—1)(-1) = 1 


Thus, Q is a rotation. 


Let u and ¢ be the unit vector parallel to the axis of rotation and the angle of rotation, respectively; 


then, 
Q = (1 — 2ee")(1 — 2ff7) = 1 — 2ee” — off? + 4(e7 f)ef? 
Hence, 
1 
vect(Q) = —54(eT fe x f =using 
and 
tr(Q) = 3 — 2 — 2 + 4(eTf)? = 4(eTf)? — 1 = 1 + 2cos¢ 
Therefore, 
sing = ||vect(Q)|| = 2le7f] lex fl, 0<ġ< 7 
whence 


exf . T 
u = ———— _sign(e’ f) 
lle x fli 


where sign(e7f) is the signum function, which is +1 if e7f > 0, —1 if eTf < 0, and undefined if 
e’f = 0. Hence, u is undefined if eff = 0, which indicates that Q = 1 — 2ee™ — 2ffT, in which case 
Q is symmetric, as in Exercise 2.13, and 


u=exf, @=7 
Hence, if e7 f Æ 0, 


Qlef f 
cos = 2(e7f)?—1, = $= tan | Jee Ee ' 


-eT 
2.9 The product Q = R: R2 is readily computed as 
Q = (1 — 2eeT)(1 — 2ff7) = 1 — 2ee? — 2ffT + 4(eTf)efT 
Let ¢ be the angle of rotation of Q, and hence, 
tr(Q) = 3 — 2 — 2 + 4(eTf)? = 4(eTf)? — 1 = 1 + 2cos¢ 


If eTf Æ 0, then 
cos ¢ = 2(e7f)? — 1, 


Because | cos ¢| < 1, we can write the above expression in the form 
—1 < 2(e"f)? -1<1 or 0<|e"f| <1. 
Since e and f are unit vectors, this condition is always respected for any e or f. However, 
vect(Q) = 4(e’ f) vect(ef”) = —2(eTf) e x f 


where we have recalled relation (2.60). The axis of rotation of Q should thus be parallel to e x f, which 
is the sole condition for Q to be factorable as Ri Ro. 
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2.10 (a) With A = CPM(a), the equation can be rewritten as 
(1+A)v=b 


which was defined as B in Example 2.3.2. In that example, it was shown that 1+ A is always invertible, 
and hence, v can be expressed as 
1 


oe ee Blas aD 


[(1 + |lal|?)1 — A + A?]b 
which, upon expansion, yields 


1 
= b+ ——,,(-Ab+ A’b 
v= b+ TF Tale! 


and, finally, 
—axb+b+(a’bja 


v= 
1+ |lall? 


(b) For v to be orthogonal to a, one must have vfa = 0, which means that, when “dotted” with a, 


the numerator of the above expression must vanish, i.e., 
—(ax b)”a +bľa + a‘ bjlal|? = 0 
— 
0 
which thus reduces to 
(a”b)(1 + lal?) =0 


As the second factor above cannot vanish, the condition sought is that a and b be mutually orthogonal. 


(c) For v to be orthogonal to b, a similar relation follows upon “dotting” the numerator of v with b, 
which leads to b 
— (a x b)’b+||b||? + (afb) = 0 
— 
0 


which holds if and only if the sum of two positive quantities vanishes. Obviously, this sum cannot 
vanish, and hence, v cannot be orthogonal to b. 


2.11 Let À be an eigenvalue of E, and u the corresponding unit eigenvector, i.e., 


Eu =u (1) 
Multiplying eq.(1) by E, we obtain 
E7u =u (2) 
and multiplying a second time by E, 
Eĉu = \°u (3) 
However, from eq.(2.51) 
E’ = —E (4) 
Substituting eq.(4) into eq.(3), we obtain 
—Eu = \°u (5) 


Moreover, from eq.(1) we have Eu = Xu, and thus eq.(5) reduces to 


(1 + A?)u = 0 (6) 
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However, ||u|| = 1 and hence, u 4 0. Therefore, we have from eq.(6) 
P(A) = A(1+ A”) =0 (7) 


Equation (7) is the characteristic equation of E, and P(A) its characteristic polynomial. The three 
roots of eq.(7) are readily obtained as 


A, =0, =j, à =-j 


The eigenvector u; corresponding to 1 is, obviously, e, for Ee = 0. 


Now, to prove that the eigenvectors associated with the two complex eigenvalues are complex and 
mutually orthogonal, assume that v and w are two 3-dimensional complex vectors. Moreover, the 
former is assumed to be the eigenvector of E associated with Az = j, the latter with A3 = —j, i.e., 


Ev = jv, Ew = -jw (8) 
The scalar product p of v times w is a complex number, given by 
p= v*w (9) 


where, if, for example, v = vr + jvc, with v, and ve denoting two 3-dimensional real vectors, then, 
the transpose conjugate of v is 
v* = vľ — jvľ (10) 


Moreover, the transpose conjugate p of p, being a scalar, is simply its complex conjugate, which equals 
the transpose-conjugate of the product on the right-hand side of eq.(10), namely, 


p= (v*w)* = wřv (11) 
Now, multiplying both sides of the second of eqs.(8) by j, we obtain 
w = jEw (12) 
Substituting w by the above expression in eq.(10), an alternative expression for p is found: 
p = Vv" (jEw) = jv" Ew (13) 
Likewise, multiplying both sides of the first of eqs.(8) by j, we obtain 
v = —jEv (14) 
which, when substituted into eq.(11), yields an alternative expression for D: 
P = w“ (—jEv) = —jw* Ev (15) 


Now, upon taking the complex-conjugate of both sides of eq.(13), one more expression for p is obtained, 
namely, 
Dp = —jw*(E’v) = jw*Ev (16) 


Notice that the two expressions for p, eqs.(15) and (16), differ only in the sign, i.e., 


=P p= -p 


SI 


which holds if and only if p = 0, thereby showing that, in fact, the two eigenvectors of E in question, 
v and w, are mutually orthogonal. In the next step, these two eigenvectors are found. The real and 
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imaginary parts of v were defined above as v, and ve, respectively. A similar definition is applied now 
to the real and imaginary parts of w, namely, w, and we, respectively. Substitution of v in terms of 
its two components, real and imaginary, into the first of eqs.(8) yields 


E(v, + jVc) = JVr — Ve (17) 
whence two real equations follow, one for the real, one for the complex part: 
Ev, = —ve, Ev. =v; (18) 
Next, upon multiplying both sides of the first of the above equations by E, one obtains 
E?v, = —Ev, = —v, 
where the second of eqs.(18) has been recalled, the foregoing equation thus becoming 
(E? +1)v, =0 


and, in light of the expression for the square of a 3 x 3 skew-symmetric matrix, eq.(2.39), the above 
equations leads to 
e(e’v,) =0 
—— 
0 


which implies that v, is normal to e. A similar reasoning leads to the same conclusion about ve, 
namely, that v. is also normal to e. Moreover, the first of eqs.(18) leads to 


Ve = —€E X Vr 


whence v, is also normal to vp. Therefore, any two orthonormal vectors lying in the plane normal to e 
are candidates for the real and complex components of v. A similar reasoning leads to the conclusion 
that any two orthonormal vectors lying in the plane normal to e are candidates for the two components 
of w. However, v and w are known to be mutually orthogonal. Upon imposing the condition v*w = 0, 
the relation below follows: 

vlw, +v] we t+ j(vi we — viw,) =0 


which readily leads to two real equations, one for the real part, one for its complex counterpart: 


T 
r 


T T; 


T 
V Wr HV We =0, vV We— vw, = 0 


The two above equations are verified under the conditions w, = ve and We = Vr, i.e., 
W = Ve + JVr 
In conclusion: any pair of orthonormal unit vectors normal to e constitute the real and imaginary 


parts of the eigenvector associated with A2 = 7. The same vectors, but with their real and imaginary 
roles exchanged, constitute the real and imaginary parts of the eigenvector associated with A3 = — j. 


2.12 From eq.(2.53), we have 
Q = e™ = f(E) 
where E is the cross-product matrix of the unit vector e with eigenvalues { 0, j, —j }—see Exercise 2.10. 


Let A, be an eigenvalue of E and yx, be an eigenvalue of Q. Then, by the Cayley-Hamilton Theorem, 


we have 
jig SOE = f(Ak) 
Therefore, 
n=l eae” pa = eI? 
10 
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2.13 Here we have 
tr(Q) = -1=1+42cos@ 


from which, 


Using eq.(2.49), we have 


Q+1 = 2ee” 
which gives 
1 2 —2 2 e? €1€2 €1€3 
3 —2 2 —2| =2 | e1e2 e2 €2€3 
2 —2 2 €1€3 €2€3 e2 
Thus, 
en 
ee am Era 33 
Choosing e1 = +V3/3, we finally obtain 
V3 1 
e=—|-l], ¢@=7 
3a 


2.15 We have, from Cayley’s Theorem, 
Q=(1-S$)(1+8)7' 
where 1 + S is invertible. Indeed, if we write vect(S) = s = ce, with o > 0 denoting the Euclidean 
norm of s, then S = oF, in which E is the cross-product matrix of the unit vector e. Hence, the 
eigenvalues of S are, from Exercise 2.10, 0, gj, and —o7. Consequently, the eigenvalues of 1+ S are 
1, 1+o7, and 1— øj, none of which is zero, thereby showing that 1+S is indeed invertible. Multiplying 
both sides of this equation by (1+ S), from the right, we obtain 


Q(1+S) =(1-S) 
i.e., 
(1+Q)S=1-Q 
and hence, 
S=(1+Q)"*(1-Q). 


For this factoring to be possible, then, matrix 1 + Q should be invertible, i.e., nonsingular. Therefore, 
the above-mentioned factoring is not possible if one of the eigenvalues of Q is —1, which happens iff 
the angle of rotation is 7. 


2.16 The given matrix Q is proper orthogonal and symmetric; therefore, Q represents a rotation with @ = 7 
(see also the solution of Exercise 2.13). As stated in the solution for Exercise 2.15, the Cayley factoring 
is impossible in this case. 


2.17 From Exercise 2.15, we have 
S=(1+Q)*(1—-Q) 


where Q can be expressed from eq.(2.54) in the form 


Q =1+sin gE + (1 — cos ¢)E? (19) 
> 14+Q=21+4sin¢E+ (1 —cos¢)E? (20) 
11 
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On the other hand, from the Cayley-Hamilton Theorem we can write 
(1 +Q)! = a1 + BE+ 7E? (21) 
where a, 8 and y are constants to be determined by imposing 


(1+ Q)(1+Q)'=1 (22) 


We can compute coefficients a, 6 and y using one of two approaches. 


First approach: using eigenvalues 


Eigenvalues of Q: {1, ef?, e~J*} (see Exercise 2.11) 
=> Eigenvalues of 1 + Q: {2, 1+e%, 1+ e739} 


= Eigenvalues of (1+ Q)7!: {5 nr =} 

Eigenvalues of E: {0, j, —j} (see Exercise 2.10) 

Eigenvalues of E?: {0, —1, —1} (squares of eigenvalues of E) 

Let X = al + JE + 7E? 

=> Eigenvalues of X: {a, a + j8 — y,a — 78 — 7} (respecting foregoing order) 
Equate eigenvalues of © with those of (1 + Q)7?: 


1 
a aks a CT (23b) 
. 1 
a al cae (23c) 
Substituting a = 1/2 from eq.(23a) into eqs.(23b &c), we obtain, after some simplifications: 
ej? — ej? —2jsing sind 
2j(2 + eib +e-i$)  2j(2+ 2cosø) 2(1 + cos ¢) 
Second approach: using relation (22) 
Using eqs.(20) and (21), we derive 
(1+ Q)7'(1+ Q) = (a1 + BE + YE”) [21 + sin gE + (1 — cos ¢)E?| 
= 2a1 + asin dE + a(1 — cos ¢)E? 
+28E + b sin dE? + B(1 — cos ¢)E* 
+2 yE? + ysin ¢E® + y(1 — cos ¢)E* 
From eq.(2.51), we have 
E? = -E (25a) 
E* = (1 — ee") = —E? (25b) 


12 
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Therefore, 


(1+Q)-1(1+Q) = 2a1 + [asin d + 28 — B(1 — cos $) — ysin ¢] E 
+ [a(1 — cos ¢) + Bsing + 2y — 7(1 — cos ¢)] E? 
Since (1 + Q)~1(1 + Q) = 1, we have 
2a=1 
(a — y)sinġ + B(1+cos¢) = 0 
(a — y)(1 — coso) + Bsing +2y = 0 


Solving these three equations for a, 8 and y, we obtain 


iL —sing 
= = = Ą— = 0 
2° f 2(1 + cos ¢)’ i 


Computation of S 


Substitution of foregoing values of a, 6 and y into eq.(21) leads to 


a d T sin @ 
(1+Q) =o" 2(1 + cos ¢) 


and 


$=(1+Q)"a-9) = (51-2 8) aa) 


_il 1-Q- sin d E sin b EQ 


2 1 + cos ġ 1 + cos ġ 
Using the representation of Q given in eq.(19), along with relations (25a), we obtain, after some 
simplifications, 
_ -sing 
~ L+cos¢ 


But, if we recall the tan-half identities: 


1-7? 2 
cos ġ = sind E T = tan ($) 


ae B 
1 +r?’ L472? 


then, the above expression is equivalent to 


S = — tan ($) 


Finally, let us take the vector of the two sides of the above equation 


apea 


which is, using the definition for the Rodriges vector p in this Exercise, p = tan(¢/2)e, the negative 
of this vector, i.e., 
vect(S) = —p 
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